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Abstract 

The common belief that the lightest supersymmetric particle (LSP) might be a neutralino, 
providing also the main Dark Matter (DM) component, calls for maximal detail in the 
study of the neutralino properties. Motivated by this, we consider the direct production 
of a single neutralino Xi a ^ a high/energy hadron collider, focusing on the Xi an d X2 cases. 
At Born level, the relevant subprocesses are qq — > gq — > XiQL,R and qq' — > XiXf] while 
at 1-loop, apart from radiative corrections to these processes, we consider also gg — > x^d, 
for which a numerical code named PLATONgluino is released. The relative importance 
of these channels turns out to be extremely model dependent. Combining these results 
with an analogous study of the direct P arr production, should help in testing the 
SUSY models and the Dark Matter assignment. 

PACS numbers: 12.15.-y, 12.15.Lk, 13.75. Cs, 14.80.Ly 
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1 Introduction 



The lightest neutralino state, x?> is often assumed to be the Lightest Supersymmetric 
Particle (LSP) pQ. As such it is also a candidate for the origin of Dark Matter (DM) [2|. 
This assumption has of course to be verified though, by analyzing the results or constraints 
reached by experiments trying to detect Dark Matter through direct or indirect methods 

However, even in the minimal MSSM version of the SUSY models, the large parameter 
space induces great uncertainties in the neutralino properties. So to check the consistency 
of the DM idea, it is essential to establish the neutralino properties through production at 
high energy hadron and lepton colliders. The first such possibility of neutralino produc- 
tion, will probably be through cascades at the CERN LHC (HUE]- But precious additional 
independent information from LHC could also be obtained by studying the smaller signals 
of the direct pair production, as well as the production in association with other 

sparticles in processes as {x1,Ql,r) or (XiXf)- When the LC collider, will finally be 

built, a wealth of additional information will become accessible [Zj. 

Studies of the pure QCD effects to these channels at LO and NLO have already 
appeared |H1 E]- A recent summary can be found in ^Hj, where the results of a NLO 
QCD computation are presented for various processes including neutralino production in 
association with a gluino, squark, slepton, chargino, or another neutralino. The overall 
conclusion of these computations is that at the LHC range, the pure QCD soft and 
collinear corrections always increase the LO cross section by an amount which, depending 
on the subprocess cm. energy and the masses of the particles involved, lies in the range 
of 10% to 40%. 

Also important at LHC though, turn out to be the leading and subleading 1-loop 
logarithmic (LL) electroweak (EW) corrections. Particularly for processes characterized 
by non-vanishing Born contributions, such effects show a largely universal structure with 
the leading ln 2 (s)-terms solely determined by the couplings of the known gauge bosons 
{W^i Zfj) to the external particles of the process; which in turn are fixed completely by 
their quantum numbers. The situation is different for the subleading single- In (s) terms 
though, which depend on the couplings and masses of all virtual particles, gauge or non- 
gauge, shaping up the underlying dynamics [TTJ H21 EH E] • Thus, depending on whether 
SUSY is "near by" with all MSSM sparticles below the TeV range, or some of the sparticles 
are very heavy, or even that the pure simple SM model stays correct till very high scales, 
will only affect the subleading single ln(s)-terms fH\ ITS} IT4*j. 

The most striking characteristic when comparing these EW corrections to the afore- 
mentioned pure QCD ones, is that they are of roughly similar magnitude, but have opposite 

sign nn n2i nsi nu 

Particularly for the subprocess qq —>■ XiXj contributing to the neutralino-pair pro- 
duction, these effects have been studied in [THj, where the calculation of the pure 1-loop 
process gg — > was also included. If the masses of the squarks of the lrst and 2nd 
family turn out to be very heavy, it might happen at LHC, that the gg XiXj* contri- 
bution is comparable to that of the LO process qq — > XiXp particularly at low invariant 
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masses where the gluon flux is very large. 

Additional information on neutralinos in a hadron collider could be obtained from the 
single neutralino production triggered by the subprocesses: 

1Q^X°i9 , 99^X°i9 , 9q^X° t QL,R , qtf XiXf , (!) 

where the indices now enumerate the neutralino and chargino respectively. The aim 
of the present paper is to study the physical consequences of these subprocesses using the 
same procedure as in Since different particles are involved in each of them, their 

combined study is sensitive to different aspects of the underlying model. 

For the first, third and fourth of the subprocesses in (|I]L this model sensitivity arises 
already at the Born level mainly caused by the (gaugino-higgsino) mixing matrices 1 
multiplying the basic gaugino and higgsino couplings. The relevant diagrams are shown 
in FigsHSl • This model sensitivity is further enhanced when including also the leading 
logarithmic part of the 1-loop corrections, calculated by following the procedure of [TTj . 
Thus, rather simple expressions for the amplitudes of these processes are reached, which 
apart from being very sensitive to the physical dynamics, should also be quite adequate 
for LHC energies and accuracies. 

Further model sensitivity is induced in the case of x°i9 production, by the contribution 
of the genuine 1-loop subprocess gg —>■ Xi9- The generic form of the relevant diagrams is 
shown in FigEJ On the basis of these, a numerical Fortran code called PLATONgluino 
is released, calculating da(gg — > Xi9)/dt for any set of real \x and MSSM soft breaking 
parameters at the electroweak scale [TTj. 

To explore the actual physical situation that might be realized within the SUSY ap- 
proach, typical MSSM benchmark models with real parameters are used ^ElEHIISI]- The 
LHC cross sections for proton proton collisions are then computed by convoluting the 
qq, gg, qg subprocess cross sections, with the corresponding quark and gluon distribu- 
tion functions taken from [21] . As in |15j . invariant mass and angular distributions are 
constructed, illustrations of which are given below. 

The results obtained in this paper should be useful for precise applications at LHC 
taking into account decay branching ratios and final state identifications. We will come 
back to this point in the conclusion. 

The organization of the paper is the following. In Section 2.1 and Appendix A.l, the 
general form of the Born amplitudes for qq — > Xi9 are given, together with the 1-loop LL 
EW and 2 SUSY QCD corrections to them, as well as the explicit Born expressions for 
the helicity amplitudes. The corresponding results for qg — > XiQ~L,R an d qq' — > X^xf are 
given in Sections 2.2 and 2.3, and Appendices A. 2 and A. 3 respectively; while in Section 
3, the 1-loop process gg — > is discussed. Finally, in Section 4 we discuss our results, 

1 The notation of JHj is used here. 

2 The SUSY-QCD corrections describe a special part of the complete QCD correction, intimately 
related to the SUSY dynamics. The fact that we consider them together with the EW corrections, rather 
than the pure QCD ones, is a matter of choice. See p2] for its exact definition. 
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and Section 5 presents the Conclusions. 



2 The processes qq j$g, gq -> x-^i?, ^ -> X<xf • 

The momenta, energies and masses in these subprocesses, as well as in gg — > y^g of 
Section 3, are defined as 

A(pi,Ei,mi) B(p 2 ,E 2 ,mj) , (2) 

where the masses of the incoming particles are neglected. Denoting by (p, 6 1 ) the final 
state cm. momentum and scattering angle, we have 



s = (p x + p 2 ) 2 = (gi + <?2^ 2 



t = (52 -P2) 2 = (Pi -<?i^ 2 



w = (52 - pi) 2 = (P2 - gi^ 2 



1 r -1 1/2 

' r - ' ' 2i rs /™ ™ \2i 1 7 



2v s + mf — m 2 A s + mh — mf 

2^1 2y/§ 



.s 



g 1= ^(l,0,0,l) , 92 = ^(1,0,0,-1), 
Pi = (.Eijpsinfl, O,pcos0) , p 2 = (E 2 , —psm9,0, —pcos9) . (3) 

The common characteristic of the subprocesses of the present Section, is that they all 
receive non-vanishing Born contributions determined by the diagrams in FigsOHl Since 
we neglect initial masses, the only needed vertices for calculating the diagrams for the 
first two processes, are those given by the neutral gaugino-quark-squark couplings 

AT{u L ) = - —J {ZSsw + 3Z»c w ) , Af{d L ) = - — J (Z»8w - 3Z 2 N t c w ) , 

ov ZswCw oy zswCw 

(4) 

A° R (u R ) = 2 -^Z»* , A™(d R ) = - ^Z»* , (5) 
and the corresponding chargino-gg^ R -ones 

Af(u L ) = —Z+ , Af{d L ) = -J-Z^ . (6) 
% % 

The notation of 10 is used for the neutralino and chargino mixing matrices, and % and j 
in (J3J EJ) and © , denote the neutralino and chargino index respectively. 
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For the third process qq' — > Xi Xj% determined by the three Born diagrams of Figs|2k,b,c 



one needs in addition the jy-chargino-neutralino couplings 

nWL _ ryN yN r>W R _ yN* 7- 1 1 7JV*v- (n\ 

2.1 The process — ► x^g to the LL 1-loop EW order. 

Writing this process in more detail as 

Qaiqi, Ai) g C2 (g 2 , A 2 ) -> 9i(P2,r 2 ) , (8) 

we denote by (ci,c 2 ) the color indices for (q, q) respectively, and by I the color index of 
g. The helicity amplitude is then written as Fx 1 x 2 ;t 1 t 2 ^ with color indices suppressed and 
Ai, A2, T±, r 2 denoting the helicities. The mass rrij in (JHJ) now describes the gluino mass. 

The Born level contributions to this amplitude arising from the two diagrams, in 
FigHk,b are 



C2C1 



QR 



C2C1 



P B(b) = g£ ^S^ (|) [-(? 2 )P^(x°)][«(^)P^(gi)] 
u — m= v z / c 2 ci 



u - m? fl V 2 



(9) 



c>ci 



where ()4I5|) have been used and g s denotes the QCD coupling. 

The explicit expressions of the Born helicity amplitudes F^ x TiT2 are given in (jA.lj) . 
To get full helicity amplitudes containing also the 1-loop LL EW and SUSY QCD contri- 
butions, the corrections in (|A.4|) and (jA.9|) should be added. The differential cross section 
is then obtained as 

dcose -TiE^Ts ^ |iWlT2 ' • (10) 

col, spins 

At asymptotic energies, much larger than all masses, both the dominant amplitudes 
(see ([A.lOp ). and the differential cross sections simplify considerably. 



3 We use the same notation as in |15II16| . 
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2.2 The process qg — > XiQl,R to the LL 1-loop EW order 

Writing this process as 

9ci(9i, Ai)^(g 2 ,// 2 ) X°(pi,Ti)g/ C2 (p 2 ) , (11) 

we denote by (ci, c 2 , Z) the color indices for (q, qj, g) respectively, while (J = L, R) 
determines the type of the produced squark of the first or second family. The helicities 
of initial quark and gluon, as well as the helicity of the final neutralino, are respectively 
described by Ai, /x 2 , Tl- Correspondingly, the polarization vector of the initial gluon is 
denoted as e 2 , and the full helicity amplitudes for the process is written as Fx im71 . As 
before, the kinematics are fixed by (0121), with rrij now describing the squark mass. 

The two relevant Born diagrams are shown in Fig|21 Since the incoming quark is 
massless, the squark specification by the index (/ = L, R), is uniquely associated with 
the quark helicity being (Ai = —1/2, + 1/2) respectively; this property remaining true 
at 1-loop LL level also. 

With the momenta and helicities defined by (fTTjl . the contributions from the diagrams 
in Figsl^b may then be written as 



F B(a) = gsAfiqj) /* 
s V 2 



I 





C2C1 



«(x?)iMrfi + rf2)y 2 «(?i)] 



^ B(b) = f^(l) Hx°)Piu(q)](2e 2 . P2 ) . (12) 

Ql 

The resulting Born helicity amplitudes appear in (jA.lll IA.12|) . while their asymptotic 
expressions are given in (jA.14|) . The universal and angular parts of the LL EW and 
SUSY corrections to these amplitudes are respectively given in (jA.16l IA~T7|) . 

After averaging over spins and colors, the cross sections are obtained from these am- 
plitudes by 

dajqg -> XiQi) _ g , p , 2 h „\ 

dcosO 30727TS ^ 1 Al ' M2 '" 1 ' ' 1 j 

col, spins 

The asymptotic expressions of the amplitudes including all LL EW and SUSY QCD 
corrections appear in (jA.18|) . 



2.3 The process qq' —> XiXj to the LL 1-loop EW order. 

The two contributing processes in this case, namely ud — > x°i + Xj an d du — > Xi + Xj> 
should give equal differential cross sections, because of the CP invariance valid for real 
soft MSSM breaking and /i parameters: 

da(ud -> x° + Xj) _ dajdu -> x° + Xj) ^ 
dcos^ dcos^ 
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We therefore concentrate on 

u(qi, Ai)J(g 2 ,A 2 ) ->• X?(pi,Ti) Xj(l>2,Tj) , (15) 

where the helicities and momenta are defined, so that (JSJ) keeps describing the kinematics 
with mj now being the chargino mass. The helicity amplitudes are denoted as Fx 1 x 2 - TiTj . 

The Born level contributions arise from the three diagrams in Figs|Hk,b,c caused re- 
spectively by the exchanges of a W + in the s-channel, a u^-squark in the t-channel, and 
a (i^-squark in the li-channel, and suitably analyzed as 

pijB _ oijB , rpijB . jjijB Qg\ 

A1A2 \TiTj A1A2 \TiTj ' A1A2 \TiTj ' A1A2 JTiTj ' \ / 

where the indices refer to the neutralino and chargino respectively. Note that, since 
we neglect quark masses, there are no R-squark exchange contributions. 

Defining the momenta and helicities as in ()15|) . and using (|7|). the contributions from 
the three diagrams in Figs|HK,b,c to the Born helicity amplitudes appear in (|A.19l IA.201 
IA.21)) respectively. 

At asymptotic energies, only F % 2 + h and F l] + ^ retain a non-vanishing Born contri- 
bution appearing in (|A.22|A.23|) ; while the associated EW universal, SUSY QCD, RG 
and angular LL corrections are shown respectively in (jA.24l IA.25|) . (|A.26j) . ()A.27|) and 
(IA.28IA.29I) . 

The spin and color averaged differential cross section is calculated from 

dajud -> x°j Xj) _ d °(du ~> X°iXj) (3 - 2 

dcos# ~" dcos# 384vrs^' AlA2 ' T1T2 ' ' 1 ' 

spins 



3 The one loop process gg — > %^ 

The momenta, helicities and color indices (ai, a 2 , 03) of the particles participating in this 
process, together with the polarization vectors of the gluons, are defined though 

0oi(?i,eiM) + 9a 2 (q2,e 2 (^2)) -> X°(pi,Ai) + <?a 3 (P2,A 2 ) . (18) 

The kinematics is defined in (jHJ), with m, denoting the neutralino mass and rrij the mass 
of the gluino. The helicity amplitude of the process denoted as F^°£?j^ X2 (0) , satisfies 

KlSxM = (-^K^lS^ - e ) > ( 19 ) 

because of Bose symmetry among the initial gluons. 

This process first appears at the 1-loop level, driven by the diagrams generically shown 
in Fig0] These consist of three types of box diagrams named (Bl, B2, B3), which are of 
exactly the same form as those met in neutralino pair production in an LC 77 collider, or 
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in the calculation of the reverse process of dark matter annihilation to photons 0] . In 
addition to them, there are three types of s-channel triangular diagrams (si, s2, s3), and 
two types of t-channel triangles. These diagrams have been calculated exactly and the 
results were used to construct the FORTRAN code named PLATONgluino which, after 
averaging over all spins and colors, calculates 

da{ " ? ^ in fb/TeV> 
at 

for any value of the subprocess cm. scattering angle 9 given in radians, and any set of 
real MSSM parameters at the electroweak scale. As with other related codes we have 
constructed, PLATONgluino may be obtained from [TTj . 



4 Results 

In this Section we discuss the LHC predictions for the direct production of a single neu- 
tralino associated with a gluino, squark or chargino, according to the four subprocesses 
presented in Sections 2 and 3. The predictions are valid for any MSSM model with real 
SUSY parameters. An exploration of the possible results has been made, using typical 
benchmark models ^HJ E3 l2"Uj . These benchmarks have also been used in other recent 
neutralino explorations, and their sole purpose is to help identifying the physical param- 
eters mainly affecting the neutralino production at LHC |U [7J [T3] . 

For the parton distribution functions inside the proton, we use the MRST2003c pack- 
age [21] at the factorization scale 

Q = Er± + E Ii (20) 

A complete summary of the relevant parton formulae and kinematics may be found in 
Appendix B of [To] . 

As observables we use the invariant mass distribution da/ds of the aforementioned 
subprocesses, and the cm. angular distribution da/d\ defined e.g. in eqs.(B.39,B.43) of 
[15j . The x~ variable is always taken to describe the particle accompanying the neutralino 
in the subprocess and is defined in terms of cm. variables 

1 --^costf* 

* = e%; = t-v* — T. > ( 21 ) 

1 + f^cosft* 

so that our treatment covers both, the LSP well as the case of a heavier \i ■ 

The transverse momentum distribution is not shown in any detail here, since it presents 
the same features as the mass distribution; a similar situation has already been noticed 
f° r Xi X°j production [T^] . Depending on the experiment of course, such distributions may 
also be useful. 
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Extensive sensitivity of the single neutralino production processes to the SUSY MSSM 
parameters, is observed. This is caused mainly by the dependence of the neutralino 
couplings on the percentage of their gaugino (Bino and Wino) or higgsino components, 
through the mixing matrices ^H]- The four processes in ([T} react differently to this 
percentage, as the first three are mainly controlled by the gaugino components 4 , whereas 
the fourth process depends both, on the gaugino and on the higgsino components. 

As we are especially interested in the structure of the LSP, supposedly the lightest 
Xi, which, depending on the benchmark, can predominantly be either Bino, or Wino, or 
higgsino, this explains the large sensitivity to the chosen benchmark. 

The single neutralino production processes are also very sensitive to the masses of the 
exchanged squarks. In the gluino case, the relative importance of the one loop process 
99 ~ > Xi9i is also strongly depending on the squark masses. For light squarks, this 
gives cross sections which are about a hundred times smaller than the ones from the qq 
process. But if the squarks of the first two generations become heavy, while those of the 
third remain rather light, it may turn out that kinematical regions exist where the 1-loop 
subprocess gg — > x°i9 is appreciable, compared to the Born-level subprocess qq — > Xi9] 
so that it cannot be ignored. Comparing with the gg — > XiXk treatment of [T^j, we 
should note that the gg — ► Xi9 process cannot be enhanced by resonance effects like those 
enhancing x°iX°k production. 

For what concerns the electroweak radiative corrections to the three Born processes, 
the computations of the leading and subleading logarithmic contributions show a negative 
effect, regularly increasing with the invariant mass, which is of the order of (10-20)% at 
the TeV range, as expected from [H]. This effect is comparable, but of opposite sign, to 
the analogous QCD correction[SJ 13 HO], and it should also be taken into account for a 
precise analysis. 

On the basis of our explorations, we present in FigsloTITUl below, illustrations for four 
different typical cases 5 : 

1) a gaugino (Bino)-type model for Xi, with light squark masses, SPSla [TH] . 
la) a Bino-type model with heavy squark masses 6 , SPSlaa, 

2) a higgsino type model for Xi, Xx > with light squark masses 7 , AD(fg9) [TT?] . 
2a) a related higgsino type model, but with heavy squark masses 8 , AD(fg9a) |T§] . 

4 Some dependence on the higgsino component appears also for the subprocess gg — ► \i9i caused by 
(t, 6)-loop contributions. It should be remembered though that the contribution of this subprocess, being 
of higher order, is generally suppressed. 

5 In all cases we have used Suspect2_3 to calculate the various masses |2"3] . 

6 SPSlaa is constructed from SPSla of JSIi by simply putting the high scale SUSY breaking soft 
sfermion masses of the first and second generations, at 5000 GeV. 

7 This model is extracted from Fig. 9 of (T^|. It is characterized by the high scale values: mi/2 = 
420 GeV, A = 420 GeV, tan /3 = 40, and m = 600 GeV for all scalar masses except m Hu = 600\/2 GeV. 
To preserve the predominantly higgsino nature of xli m t = 174 GeV should be used here, as this was 
the case when the model was constructed. 

8 It is constructed from the EW scale masses of AD(fg9), by only changing the sfermion soft SUSY 
breaking masses of the first two generations, which are now put at 5000 GeV. 
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discussed in turn below: 

• Let us first consider the Xi gaugino-type model (SPSla). This model gives invariant 
mass distributions for the three Born processes in (|T]). which are largely observable 
in the 1 TeV range; i.e. cross sections of about lOOfb for the first two cases, but only 
lOfb or less for the third one. It also gives an invariant mass distribution for the 
qg — > XiQL,R channel, which is more important at low masses, due to the behavior 
of the gluon distribution function. The angular distribution for qq — > x°9i described 
at Born level by the t- and u- channel squark exchanges indicated in Fig^ flattens 
out at large Xg m this model. 

It is amusing to remark that a very similar behavior is also expected in the universal 
m-SUGRA type model which has been identified by |2lj, as "a best" description of 
all present particle and cosmological constraints 

• Comparing SPSlaa to SPSla, one notices a reduction of the invariant mass distribu- 
tion in the 1 TeV range, by more than an order of magnitude for by somewhat 
less than an order of magnitude for XiXii an d, obviously, a complete suppression of 
the x®*! production. Moving from SPSla to SPSlaa, there appears also a change in 
the ^-distribution, which becomes steeper for but remains roughly similar for 

• We next turn to the higgsino type model AD(fg9) for the three channels studied 
here. Comparing them to the gaugino SPSla and SPSlaa models, we find that the 
predicted cross sections for the x®9 an d XuQ channels, are much smaller now, since 
they are controlled essentially by the gaugino components; (see Figs l8ll(jj) . On the 
contrary, the cross section for the XiXi process is larger, because of the presence 
of the s-channel W exchange diagram involving higgsino components; (see FigOJ). 
The ^-distribution becomes also flatter, for the same reason. These can be seen by 
comparing FigslHtEt andlTUb.d with Figs^tEfc and Ho, d respectively. 

• Finally, we compare the results of the models AD(fg9) and AD(fg9a), in both of 
which there is a large higgsino component to x\i as wen as t° X2- 111 AD(fg9a) the 
squark contribution, coupled through the gaugino component, is further reduced 
compared to AD(fg9), leading to an even smaller prediction for the Xi9'i compare 
FigsJEHEH with fTOk.c. On the other hand, for the X1X1 an d X2X1 channels, the 
cross sections are comparable, since they both receive a large contribution from the 
higgsino component coupled to the intermediate W boson, which is not affected by 
the change in the squark mass; compare FigslHtEt withfTUb.d. 

Finally we comment on the difference between the magnitude of the x\ production 
cross section (in which we are mainly interested in), and the x® one ! X2 production is 
generally more copious than the Xi one, becoming progressively more pronounced as we 
go from the Xi9 channel, to and eventually to the x'iXi channel, where it reaches a 
factor 10 in the SPSla model. This factor is even larger in the SPSlaa model; of order 
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100. These differences are due to the mixing matrix elements appearing in the Born 
amplitudes, which control the Bino, Wino and higgsino components of the neutralinos. 



5 General Conclusion on \i production 

In this paper we have analyzed the single neutralino production processes Xi9i Xi9 an d 

x°xt at LHC - 

The complete set of helicity amplitudes for the Born terms of the subprocesses qq — > 
Xi9i 91 ~* XiQL,R, qq' —> XiXj nas been written down, together with the leading and 
subleading logarithmic electroweak corrections to them. Compact analytic expressions 
are presented, which are applicable to any MSSM model with real parameters. We have 
also included the complete one loop calculation of the subprocess gg — > xWj f° r which a 
numerical code called PLATONgluino is released [T7j- 

The pure QCD corrections, which have already been given in previous papers |§1 110j. 
have not been reexamined. But we have emphasized that, contrary may be to naive 
expectations, the leading logarithmic EW and QCD corrections at the LHC energies have 
similar magnitudes, but opposite sign. Thus, they should both be taken into account in 
analyzing the experimental data. 

The single neutralino production processes have been found to be mainly sensitive on 
two physical sets of quantities; namely the amount of gaugino and higgsino components of 
the neutralinos, and the scale of the soft breaking parameters for the squarks of the first 
and second generations. To emphasize this, a set of illustrations for LHC invariant mass 
and angular distributions have been presented which indeed show this sensitivity. These 
were based on four "benchmark" models, but more were explored in our actual runnings. 

This physics output should of course be joint to the one that can be obtained from 
the XiX^ production studied previously [TS]. For that purpose we have added Figs ll 11121 
which show the invariant mass and x distributions for X2X1 an d X2X2 production, in the 
same SPSla and the SPSlaa models used for the single neutralino case 9 . In going from 
SPSla to SPSlaa, one sees a reduction of the Born contribution, rather similar to what 
happens in the Xi 29 case; but one also sees that the relative role of the one loop gg process 
in SPSlaa, is more important for X2X1 production, then for Xi29- So the neutralino pair 
production channel has its own typical features. 

Summarizing, we have observed that the channels xl,9i Xi9i XiXj an d XiXj present an 
important sensitivity to the neutralino structure; particularly to the relative magnitude 
of its gaugino and higgsino components. They also present a considerable sensitivity to 
the MSSM mass spectrum for the gluino, squarks, charginos and Higgses; the later being 
able to lead to possible resonance effects. 

We conclude by emphasizing that the results obtained in this paper should be com- 
pleted by detail experimental studies dedicated for LHC. Observables should then be 

9 Analogous results have been shown in Figs. 11, 12 of ^H]j where though, \ is defined as the inverse of 
the present one. 
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constructed addressing neutralino, gluino and squark decay channels to various numbers 
of jets and leptons. Such observables should also reflect, at some important level, the 
sensitivity to the neutralino properties 10 , that we have observed at the level of the basic 
processes. 

We hope that their measurement will be able to confirm or infirm the possibility that 
the neutralino is an important component of the Dark Matter of the Universe. 
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10 If discovered at LHC or a Linear Collider, it would also be the first time that physical quantities 
highly sensitive to the Majorana nature of a particle reach such a high observability. 
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Appendix A.l Helicity amplitudes for qq — > y®g 

Starting from the explicit expressions of the Born helicity amplitudes for the process 
shown in ©, are 
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+ 
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+ 
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Cg sin 9 
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C g sin 6 
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+ + TT 1 



(A. 



where 



9s (X 1 



(A.2) 



C2C1 



with (ci, c 2 , I) denoting the color indices of the quark, antiquark and gluino respectively, 
as defined in (jSJ). The kinematics is determined in 

In the high energy limit, where s, \t\, \u\ are all much larger than all masses, the only 
non-vanishing Born amplitudes simplify to 



j-iB asym 



t-,B asym 
+ 

pB asym 

FB asym 
j l'_ 



-g s V2Af{q L 



g s V2Ar(q R )(j 



C2C\ 



7 N T 



(3) 



= 2eg s 

L6c^ sw 

2eg B ®±Z»*(f 
cw V2> 



9 rvN 



C2C\ 



-1eg s 



C2C1 

I 



ZN* T 
li , J 9 



(3) 



+ ■ 



-g s V2Ar (&) (|) = -2e^S (f ) 

V Z / C2C1 C^/ \ Z / c 2 ci 



7N* 
J 2i 



(A.3) 



The 1-loop universal leading logarithmic EW and SUSY QCD corrections only affect 
these asymptotically dominant amplitudes and are given by [TTJ [T2J EH E] 

n- gL + {2iiy-{wfM^) z» , 



zpUniv 


j-iB asym 
— *-++- 


iriUniv 
*-+-+ 


= F B 


asym 
h- + 


TjiUniv 
*+—+ 


= *?- 


asym 


zpUniv 


= 


asym 
-+- 



V 2 



C2C1 



•7? 7;L -(2/J)c-(W)^f^) Z 2l . 



C2C1 



K- gR , (A.4) 

where 21^ = ±1, depending on whether q = u or d. The universal LL correction due to 
the gg-pair is contained in the parameter 

K~ gH = c ew (qq, gauge)^ + c SQCD (qq) (A.5) 
in (I A. 41) where H = L or R, and 



s 



describes the SUSY QCD correction, while 



a 



r'"( W . W u W 0// = _^ J ^ + 1) + j? ( 2hj _l__ hj ^__) . (A.7) 
Air I s A w 4c^J - - 



S , 2 S 
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gives the purely EW one. Here (I q , Y q ) = (1/2, 2/6) should be used for H = L; and 
(Jq, Yg) = (0, 2Q q ) for H = R, with Q q being the quark charge. Since we neglect 
quark masses, the associated Yukawa contributions are also neglected in (|A.4J) . which is 
legitimate for the quarks found as partons inside the proton. 

Finally, the universal correction due to the final neutralino appearing in (jA.4|) , is given 
by [HI EE2 EH 



'{W) = c ew (W) 



a 



^s 2 w 



M 2 



(A.8) 



which is solely induced by the Wino component of the neutralino. 



The only other EW correction that appear within the 1-loop LL level, is the angular 
one, given by PH Q2J ECU 



pang 
^ang 
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rn 
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s s . 

, ~t , — & 
In h In 

s s . 



'1 

2 / C2C1 

t) 

Z / C2C1 



(A.9) 



which arises solely from gauge exchanges between the neutralino line, (of which only the 
Wino component contributes), and either the quark or the antiquark lines [TT| I12[ ITS]. 
No one loop EW Renormalization Group (RG) corrections are generated in this case 

By adding to the Born helicity amplitudes in (|A.1J) . the corrections (|A.4|) and (|A.9J) . 
the complete helicity amplitudes are constructed, including all LL 1-loop electroweak 
effects. 

Taking into account all above corrections at asymptotic energies, the dominant am- 
plitudes may be written as 
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(^) 2eg s ^Z»[l + K SR ] . 



(A.10) 
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Appendix A. 2 Helicity amplitudes for qg — > XiQL,R 

Defining momenta and helicities as in the Born amplitudes in (II 2J) lead to the helicity 
amplitudes 



F B 



(q c i9l -> X%c 2 ) = FfSr^ciQl -> x"<?/c 2 ) + F^2;rM<*9l -> X??/ 



,B(b) 



C2^ J 



(A.11) 



with the separate contributions from the two diagrams in Figs|2k., D giving 

p B(a) ____ ^ 



Ai/i2;n 



2s 1 /4 v /£ 



==(4-) {(^i + +P)[^2(Ai - n) + (Ai - n) 2 ] cos 



+ (E 1 + m, L -p)[(Ai +Ti) +/i 2 (Ai + r 1 ) 2 ]sin-|5/ Al 



^B(b) 
AiAi2;ri 



</.4 9/ (?/) s 3 / 4 /3/i 2 sin^/A 



2 V /(E 1 + (t-mf,) 
9 ~) 

{E l +mi + p)(Ai - ri) Edn-jtf/Ai 



(— ) {( J B 1 + m i -p)(A 1 + n) 2 cos 

\ Z / C2C1 v 



(A.12) 



where by a slight abuse of notation, 5i\ 1 simply indicates that (/ = L, R), is uniquely 
associated with the quark helicity being (Ai = —1/2, + 1/2) respectively. An alternative 
expression might be obtained by substituting in (|A.12|) 



Ei + TUi ± p 



mi 



s + mi) 2 — m 2 ± 



s — rrii) 2 — in 2 - 



(A.13) 



At asymptotic energies (much higher than all masses), (|A.11| IX.12J) imply that there 
is only one non- vanishing amplitude for each squark type; i.e. 



-iB asym/ 



"l9l —> XiQLc 2/ 



-g s V2A° L (q L 



A' 



2eg s 



7 f 



N 



+ 



Qc w s w 



r(3) 



COS 

C2ci V2 

AS (9 
cos (- 

C2Ci \Z 



-^^A 0il fe)§) c2ci cos(| 



(A.14) 



with the neutralino squark couplings in the r.h.s. determined by (JUEJ). 

As in the case of |Appendix A.l| the 1-loop Universal EW and SUSY QCD LL cor- 
rections, only affect the dominant amplitudes in (|A.14[) . and they are associated to the 
quark, squark, or Wino component of the final neutralino line. Defining now 



K ¥ = c cw (qq, gauge): +-[c s Q CD (qq) + c s « CD (q|)] 



(A.15) 
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where (|A.7IA.6|) are used and the Yukawa terms have again been neglected, the net 
resulting universal correction is 



F^l-Mg^m) = - mi)** 



+ c™(w)(2i:)z»r^)r-) cos d)s^ 



s w 



C2C1 



(A.16) 



where (I A. 71 IA.81 IA.14D are used for I = L or R. 

The one loop angular electroweak corrections are induced by inserting to the diagrams 
in Fig|21 either a ly-exchange between the neutralino and the g-leg, or a (W, B) exchange 
between the q and q legs. The first case induces a ln 2 (— t)-term, while the second a 
ln 2 (— u)-one. The net result is 



Fl^ + (q cl gi -> XiQL 
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eag s /A' 
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[ln 2 (-u) -ln 2 (s)] 
i CosQz i t[ln 2 (-n)-ln 2 (5)] , (A.17) 



and as in |Appendix A.l[ there are no one loop EW RG corrections. 

By adding to the Born helicity amplitudes in (jA.llJ) . the corrections (|A.16J) and (|A.17|) . 
the complete helicity amplitudes are constructed, including all LL 1-loop electroweak 
effects. At asymptotic energies these acquire the form 
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(A.18) 
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Appendix A. 3 Helicity amplitudes for ud —> XiXj 



For the process in ()15|1 . the respective contribution to the helicity amplitudes in 1)16)1 from 
the three diagrams in FigsOk,b,c are 
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(A.19) 
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(A.20) 
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ijB 
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(A.21) 



At energies much higher than all masses, only two non vanishing helicity amplitudes 
remain, which simplify to 



pijB, as _ gijB, as _|_ jjijB, as 



rpijB, as 
r -++- 



ijB, as j_ rpijB, as 



Qll-JD, as _|_ rp 



(A.22) 



with 



s 



-++- — 



y/2i 



W 

;i - cos9)0^ L = - cos 9)[z^zt;-^=z^* 



w 



V2i 



w 



V2 



T jb^ „ a ol { ~ l)A cl* { ~ l) = _^_( Z n 8w + 3z? iCw )zt; , 



(Z^*s w - 3Z£*c w )Z lj 



N*, 



As before, the Universal 1-loop, purely gauge, EW LL corrections are 
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The first term in both expressions (jA.24l IA.25|) are due to the quark external lines, the 
second and third come from the s-channel diagram in FigJ2^, while the last term in (|A.24|) 
and (jA.25|) are induced from the u- and t-channel diagram in FigsOfc and b, respectively. 

Finally, the third term in both (jA.24|) and (jA.25|) . is a Yukawa contribution induced 
by the higgsino components of the chargino and neutralino produced through the diagram 
in FigOt. This Yukawa contribution appears in (jA.24IA.25|) . in spite of the fact that the 
participating quarks are massless jllj . 

The SUSY QCD 1-loop universal LL corrections are 
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-iijB as 
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(A.26) 



while, in this case, there exist also a 1-loop RG single-log contribution, caused by the 
W^ + -exchange in Fig|2K, which is ^T] 
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Finally the 1-loop LL EW angular corrections are 
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(A.27) 



(A.28) 



(A.29) 



The fully corrected helicity amplitudes are obtained by adding to the Born expressions 
$KM\EM\EM , the corrections $EM\EM, (EH , (!A~27j) and (|A~28l lA~29l) . 
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Figure 5: SPSla s-distributions for X12 production, in association with either g or q L R 
or xt; {Xj = xf)- 
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Figure 6: SPSla ^-distributions for X12 production, in association with either g or q L R 
or xt; {Xj = xf)- 
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Figure 10: AD(fg9a) s- and ^-distributions for Xi,2 production in association with either 
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Figure 11: SPSla s- and ^-distributions in XiXj production for j — 1 and j = 2. 
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Figure 12: SPSlaa s- and ^-distributions in XiXj production for j = 1 and j = 2. 
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